Quasiperiodic graphs: structural design, scaling and entropic properties 
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A novel class of graphs, here named quasiperiodic, are constructed via application of the Horizontal 
Visibility algorithm to the time series generated along the quasiperiodic route to chaos. We show how 
the hierarchy of mode-locked regions represented by the Farey tree is inherited by their associated 
graphs. We are able to establish, via Renormalization Group (RG) theory, the architecture of 
the quasiperiodic graphs produced by irrational winding numbers with pure periodic continued 
fraction. And finally, we demonstrate that the RG fixed-point degree distributions are recovered via 
optimization of a suitably defined graph entropy. 
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Quasiperiodicity is observed along time evolution in 
nonlinear dynamical systems [IHjI and also in the spa- 
tial arrangements of crystals with forbidden symmetries 
d, Q- These two manifestations of quasiperiodicity 
are rooted in self-similarity and are seen to be related 
through analogies between incommensurate quantities in 
time and spatial domains 5]. Here we point out that 
quasiperiodicity can be visualized in a third way: in the 
graphs generated when the Horizontal Visibility (HV) al- 
gorithm [g, 0] is applied to the stationary trajectories of 
the universality class of low-dimensional nonlinear iter- 
ated maps with a cubic inflexion point, as represented by 
the circle map Q. 

The idea of mapping time series into graphs has been 
presented in recent works |8l-[l4| where different ap- 
proaches have been developed. In particular, the period- 
doubling bifurcation cascade has been analyzed in the 
light of the HV formalism [1^ [ig and a complete set of 
graphs, called Feigenbaum graphs, that encode the dy- 
namics of all stationary trajectories of unimodal maps 
has been provided. The Feigenbaum scenario is one 
of the three well-known routes to reach chaos in low- 
dimensional dissipative systems (along with the inter- 
mittency route and the quasiperiodicity route) [l|-|3|. In 
this Letter we characterize the structural, scaling and 
entropic properties of the graphs obtained when the HV 
formalism is applied to the quasiperiodic routes to chaos. 
As we shall see, a Renormalization Group (RG) treat- 
ment of such graphs is the instrument that grants access 
to our main results. 

We briefly recall that the standard circle map [1-3] is 
the one-dimensional iterated map given by: 

0t+i = h.KiOt) ^0t+n-^ sm{27:9t), mod 1, (1) 

ZTT 

representative of the general class of nonlinear circle 
maps: 0t+i = fn,K{Ot) = dt + ^ + K-g{6t), mod 1, where 
g{6) is a periodic function that fulfills g{9 + 1) = g{9). 
The HV graphs obtained for this family of maps exhibit 
universal properties that without loss of generality we 
explain in the next paragraphs in terms of the standard 



6n+i = en+n; mod 1 
CO = n = 5/8 

2/8 




n - 














^ 


T 








T 










- 




































'■ 


































T 













^^S^^^ 



8 






1 






8 

8 




1 > 






8 

a 


^ 









[0 = 5/8 n = 0.61 17365289... 



FIG. 1: Examples of two standard circle map periodic series 
with dressed winding number u) = 5/8, 7^ = (top) and 
K = 1 (bottom). As can be observed, the order of visits on 
the circle and the relative values of On remain invariant and 
the associated HV graph is therefore the same in both cases. 



circle map. 

The dynamical variable < 6*4 < 1 can be inter- 
preted as a measure of the angle that specifies the tra- 
jectory on the unit circle, the control parameter 51 is the 
so-called hare winding number, and K in a, measure of 
the strength of the nonlinearity. The dressed winding 
number for the map is defined as the limit of the ra- 
tio: w = limt_).oo (6*4 — 9Q)/t and represents an averaged 
increment of 9t per iteration. For Q < K < 1 trajecto- 
ries are periodic (locked motion) when the correspond- 
ing dressed winding number a;(i7) is a rational number 
p/ q and quasiperiodic when it is irrational. The wind- 
ing numbers uj{^) form a devil's staircase which makes a 
step at each rational number uj = p/q and remains con- 
stant for a range of 57. For K = \ {critical circle map) 
locked motion covers the entire interval of 51 leaving only 
a multifractal subset of 51 unlocked. 

The resulting hierarchy of mode- locking steps a.t K = 1 
can be conveniently represented by a Farey tree which or- 
ders all the irreducible rational numbers p/q G [0,1] ac- 
cording to their increasing denominators q. In the devil's 
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FIG. 2: Six levels of the Farey tree and the periodic motifs of the graphs associated with the corresponding rational fractions 
pjq taken as dressed winding numbers cj in the circle map (for space reasons only two of these are shown at the sixth level) . (a) 
In order to show how graph concatenation works, we have highlighted an example using different grey tones on the left side: 
as 1/3 > 1/4, G(l/3) is placed on the left, G'(l/4) on the right and their extremes are connected to an additional link closing 
the motif G(2/7). (b) Five steps in the Golden ratio route, b = 1 (thick solid line); (c) Three steps in the Silver ratio route, 
h — 2 (thick dashed line) . 



staircase, w(r2), the width of the steps (intervals where 
w is constant) becomes smaller when the denominator q 
increases. Furthermore, if we have two steps with wind- 
ing numbers "pjq and p' /q' ■, the largest step between them 
has a winding number (p+p')l(q + q'), which is also the 
irreducible rational number with the smallest denomina- 
tor. Thus, the Farey tree also orders all mode-locking 
steps with Lo = p/q in the circle map according to their 
decreasing widths 17]. 

The HV algorithm assigns each datum 9i of a time 
series {^i}i=i,2,... to a node i in its associated HV graph, 
and i and j are two connected nodes if 6i , Oj > On for all n 
such that i < n < j. Without loss of generality, we apply 
the HV algorithm to the superstable orbits of the critical 
circle map {K = 1) with an irreducible rational number 
u}(fl) = p/q. Thus, the associated time series always 
contains ^o = as one of its values and has period q (cf. 
[ll-Q). The associated HV graph is a periodic repetition 
of a motif with q nodes, p of which have connectivity 
k = 2. (Observe that p in the map indicates the number 
of turns in the circle to complete a period). For K < 1, 
the order of visits of positions in the attractors and their 
relative values remain invariant for a locked region with 
(jj = p/q 113], such that the HV graphs associated with 
them are the same. In fig. [1] we present an example 
where the first and the last node in the motif correspond 
to the largest value in the attractor. 

In fig. [5] we depict the associated HV periodic mo- 
tifs for each p/q in the Farey tree. We observe straight- 
forwardly that the graphs can be constructed by means 
of the following infiation process: let p/q be a Farey 
fraction with 'parents' p'/q' < p"/q", i.e., p/q = {p' + 
p")/W + <?")■ The 'offspring' graph G{p/q) associated 
with w = p/q, can be constructed by the concatenation 
G{p" /q") © G{p'/q') of the graphs of its parents. By 
means of this recursive construction we can systemat- 



ically explore the structure of every graph along a se- 
quence of periodic attractors leading to quasiperiodic- 
ity. A standard procedure to study the quasiperiodic 
route to chaos is fixing K — 1 and selecting an irra- 
tional number uioo S [Ojl]- Then, a sequence a;„ of 
rational numbers approaching ujoo is taken. This se- 
quence can be obtained through successive truncations 
of the continued fraction expansion of Woo- The corre- 
sponding bare winding numbers ri(aj„) provide attractors 
whose periods grow towards the onset of chaos, where 
the period of the attractor must be infinite. A well- 
studied case is the sequence of rational approximations 
of LJoo = 0"^ = (\/5 - l)/2 ~ 0.6180..., the recipro- 
cal of the Golden ratio, that yields winding numbers 
{uJn = Fn-i/Fn}n=i,2,3... where Fn IS the Fibonacci num- 
ber generated by the recurrence Fn = Fn-i + Fn-2 with 
Fo = 1 and Fi = 1. The first few steps of this route 
are shown in fig. Hfb): cji = 1/1, 002 = 1/2, uj^ — 
2/3, a;4 = 3/5, wg = 5/8..., wg = 8/13... . Within the 
range il{Fn~i / Fn) one observes trajectories of period F„ 
and, therefore, this route to chaos consists of an infinite 
family of periodic orbits with increasing periods of val- 
ues Fn, n — )■ cx). If we denote by G^-i{n) the graph 
associated to a;„ = Fn-i/Fn in the Golden ratio route, 
it is easy to prove that the associated connectivity dis- 
tribution P{k) for G^-i{n) with n > 3 and fc < '^ + 1 
is P„(2) = Fn-2/Fn, P„(3) = Fn-3/Fn, P„(4) = and 
Pfiik) = Fn-k+i/Fn- In the limit n — >■ cx) the connectiv- 
ity distribution at the accumulation point G^-i{oo), the 
quasiperiodic graph at the onset of chaos, takes the form 



(2) 



Figure[3]^a) shows that the theoretical degree distribution 




of the quasiperiodic graph for the route described above 
is in perfect agreement with the same quantity obtained 
by applying the HV algorithm to a circle map time series 
with a dressed winding number Woo = 4'~^- The pro- 
cedure explained for the Golden ratio can be repeated 
for the 'time reversed sequence': {w„ = Fn-2/Fn = 
1~ Fn-i/ Fn}n=i^2,3...- lu this casc the ratio converges to 
1 — 0^^ in the limit n — > 00. The connectivity distribu- 
tions of the graphs {Gi_0-i(n)}„=i 2,3,.. are the same as 
in the Golden ratio route because these graphs are sym- 
metric mirror versions of the former (we use the term 
'time' because the 'time reverse' of a graph from a series 
generated by a clockwise rotation in the circle map corre- 
sponds to the graph from the same but counterclockwise 
rotation). 

The previous results can be interpreted through a 
suitably-defined Renormalization Group (RG) transfor- 
mation. We proceed as in previous works [la, [la] 
and define the RG graph operation TZ as the coarse- 
graining of every couple of adjacent nodes where one of 
them has degree k = 2 into a block node that inher- 
its the links of the previous two nodes. If we continue 
with the case of the Golden ratio, we first note that 
n{G^-i{n)} = Gi_^-i(n- 1) and TZ{G,_^-i{n)} = 
G^-i{n — 1), so the RG flow alternates between the 
two mirror routes. If we define the operator 'time re- 
verse' by G^-i{n) = Gi„0-i(n), the transformation be- 
comes TZ{G^-i{n)} = G^-i{n—l) and 'R-{Gi_^-i{n)} — 
Gi_0-i(n— 1). Repeated application of 7^ yields two RG 
flows that converge, for n finite, to the trivial fixed point 
Go (a graph with P(2) = 1). The accumulation points 
n — >■ cx), the quasiperiodic graphs, act as nontrivial fixed 
points of the RG flow: 7^{G0-i(oo)} — G^-i{oo) and 
7?.{Gi_0-i(oo)} = Gi_0-i(oo). 

The above RG procedure works only in the case of the 
Golden ratio route. (As a counterexample look at the so- 
called Silver ratio route shown in fig. ^). To extend the 
above formalism to other irrational numbers, we develop 
the following explicit algebraic version of TZ and apply it 
to the Farey fractions associated with the graphs. 
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along with the algebraic analog of the 'time reverse' op- 
erator R{x) = 1 — R{x). Observe that along the Golden 
ratio route fractions are always greater than 1/2 and we 
can therefore renormalize this route by making 
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x is x^ -|- X — 1 = 0, 



whose fixed-point equation R(x) 
with 0~^ a solution of it. 

A straightforward generalization of this scheme 
is obtained by considering the routes {aj„ — 

Fn-l/F.n}n=l,2,3... with F„ = hF„_^ + F„_2, Fq = I, 

-Fi = 1 and b a natural number. It is easy to see that 



\\mn^aoFn-i/Fn = {—b + y/P + 4:)/2, which is a solu- 
tion of the equation x^ + bx — I = 0. Interestingly, all 
the positive solutions of the above family of quadratic 
equations happen to be positive quadratic irrationals in 
[0, 1] with pure periodic continued fraction representa- 
tion: (j)^^ — [b,b,b,...] — [b] {b = 1 corresponds to 
the Golden route). Every b > 1 fulfills the condition 
F„_i/F„ < 1/2, and, as a result, we have 
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= R^ 



Fn-l 



F„. 



[b - 1)F„ + F„-2 



(5) 



The transformation Ri can only be applied (6— 1) times 
before the result turns greater than 1/2, so the subse- 
quent application of R followed by reversion yields 
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It is easy to demonstrate by induction that 
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whose fixed-point equation R('^^{x) = R2[Ri (a;)] 
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leads in turn to : 



bx—l — 0, with (j}f^ a solution of it. 



We can proceed in an analogous way for the symmetric 
case aj„ = 1 — (F„_i/F„) but, as the sense of the inequal- 
ities for 1/2 is reversed, the role of the operators Ri and 
R2 must be exchanged. 

The previous result indicates that graphs must be 
renormalized via TZ''{Gs-i{n)} — G,-i(n— 1). Again, 
the iteration of this process yields two RG flows that 
converge to the trivial fixed point Go for n finite. The 
quasiperiodic graphs, reached as accumulation points 
(n — ;> cxi), act as nontrivial fixed points of the RG flow 
since^{G^-i(oo)} = G^-i(oo). 

We observe that for fixed b > 2, and from the con- 
struction process illustrated in fig. [2ja), it can be de- 
duced that Foc(2) = 0^\ Foo(3) = I - 20"^ and 
P^{k ^ bn + 3) ^ Oyn e N. Poo{k = bn + 3), n e N 
can be obtained from the condition of RG fixed-point in- 
varance of the distribution, as it implies a balance equa- 
tion Poo(fc) = (f'b^Pooik + b) whose solution has the form 
of an exponential tail. The degree distribution Poo(fc) for 
this quasiperiodic graphs is therefore 



Poo(fc) = 



I - 20, 



k = 2 



{I - (t)b^)4>'i^ ^^'^ fc = 6n-h3, neN 



(8) 







otherwise. 



A perfect agreement between theoretical and numerical 
results for some examples can be observed in fig. Eljb). 

Notably, all the RG flow directions and fixed points 
described above can be derived directly from the in- 
formation contained in the degree distribution via 
optimization of the graph entropy functional H = 
— ^^2 ^(^) logP(fc)- The optimization is for a fixed 
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FIG. 3: Empty circles stand for theoretical degree distribu- 
tions of quasiperiodic graphs. Filled values have been ob- 
tained by direct application of the HV algorithm to critical 
circle map series of lO'' values. Distributions have been shifted 
from each other to enhance visualization. From down-up: (a) 
cjoo = [1], n = 0.606661.... (b) ojoo = [2], Q. = 0.418864...; 
cjoo = [3], n = 0.323873...; tJoo = [4], fi = 0.271502.... 



b and takes into account the constrains: P(2) — (pj^^, 
P{3) = 1 — 2(f>^ , maximum possible mean connectivity 
(/c) = 4 [H] and P{k) = 0Wk^bn + 3, n £ N. The de- 
gree distributions P{k) that maximize H can be proven 
to be exactly the connectivity distributions of equations 
[5] and m for the quasiperiodic graphs at the accumula- 
tion points found above. This establishes a functional 
relation betvireen the fixed points of the RG fl.ovi' and the 
extrema of H as it was verified for the period-doubling 



route [il,[3. 

We have demonstrated the capability of the HV algo- 
rithm for transforming into graph language the univer- 
sal properties of the route to chaos via quasiperiodicity 
in low-dimensional nonlinear dynamical systems. The 
outcome is a novel type of graph architecture where the 
motifs are the building blocks with which quasiperiod- 
icity is expressed recursively via concatenation. Signifi- 
cantly, the HV formalism leads to analytical expressions 
for the degree distribution, a function that in all mode- 
locking regions is essentially exponential. The networks' 
scaling properties can be formulated in terms of an ad 
hoc RG transformation for which the nontrivial graph 
fixed points capture the features of the quasiperiodic ac- 
cumulation points. As we have seen, it is through the 
properties of the RG transformation presented above that 
the relevant details of the quasiperiodic graphs studied 
are determined. This class represents all the quasiperi- 
odic attractors reached when irrational winding numbers 
with pure periodic continued fractions are used as dressed 
winding numbers. Furthermore, a graph entropy is in- 
troduced via the degree distribution and its optimization 
reproduces the RG fixed points. 

By means of the HV algorithm, we have found a con- 
nection between pure periodic continued fractions and 
the degree distribution of their associated quasiperiodic 
graphs. It seems feasible to generalize our results beyond 
to periodic continued fractions or any irrational with a 
pattern in its continued fraction. Finally, it has not es- 
caped our notice that, as we have a one-to-one correspon- 
dence between graphs and rational numbers, a possible 
graph algebra can be explored. 
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